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We describe a simple procedure to estimate the high-energy neutrino flux from the observed 7-ray spectra 
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of galactic cosmic ray sources that are transparent to their gamma radiation. We evaluate in this way 
the neutrino flux from the supernova remnant RX J1713. 7-3946, whose very high-energy 7-ray spectrum 
(assumed to be of hadronic origin) is not a power law distribution according to H.E.S.S. observations. The 
' corresponding muon signal in neutrino telescopes is found to be about 5 events per km 2 per year in an 

£SJ ■ ideal detector. 

PACS: Neutrinos from CR 98.70. Sa; Observations of -/-rays 95.85. Pw; SNR in Milky Way 98.38. Mz. 
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r*""- ■ The recent observations of 7-rays above TeV by H.E.S.S. are of great interest [HE]- They 

will certainly help in answering the old question of the origin of the cosmic rays till the 
knee [31 IU 03 El U\ and at the same time they could provide us a reliable guidance for what 

"pS we should expect in neutrino telescopes, at least for certain sources. 

Qh' This is evident for the main candidate sources of galactic cosmic rays, supernova rem- 

nants (SNR) [HI OH- The huge kinetic energy of the gas of the SNR could be effectively 
converted into cosmic rays by diffusive shock acceleration |lf)| . producing enough cosmic 
rays to compensate the losses from the Milky Way. When the SNR is surrounded by 
matter that can act as a target for cosmic rays, we would have a point source of very 
high-energy (VHE) gamma radiation, which seems in agreement with certain observa- 
tions. Since we expect that the matter around SNR is not too dense anyway, the 7-rays 
are not significantly absorbed, and there is a rather direct relation between VHE 7-rays 
and neutrinos. 

More in general, we think that it is important to take advantage of the detailed ob- 
servations of 7-rays whenever they exist in order to formulate definite expectations for 
neutrino telescopes. This is certainly true after the most recent H.E.S.S. observations, 
that are beginning to find VHE 7-ray spectra that deviate from power law distributions 
above 10 TeV or so (2]. 

Our recipe to calculate the neutrino fluxes is described in Sect. and the application 
to the SNR RX J1713. 7-3946 is in Sect. [31 In essence, these results are a straightforward 
application of standard techniques (and we follow as much as possible the conventions 
of [I] to emphasise this fact) but we hope that they are useful in the present moment, when 
the high-energy gamma astronomy is flourishing and the neutrino telescopes are finally 
becoming a reality. 



1 



2 Deriving the neutrino flux from the 7-ray flux 



Let us assume that the VHE 7-ray flux Fj observed from a certain source is of hadronic 
origin and that it is not significantly absorbed-the source is 7-transparent. 1 From the 
well-known relation F-y(E) = JJ° dE' 2 F w o(E')/E' valid for high-energy 7-rays we find: 

= -f § (i) 

that implies that the 7-ray flux has to be strictly decreasing. This equation, together with 
the approximate isospin-invariant distribution of pions: 

F w = F n o « F n+ « F n - (2) 

permits us to predict the flux of neutrinos using the observed 7-ray flux. It is important 
to note that the charge asymmetry has a small or negligible impact on the observable f M 
flux, compare ^JE]- The flux from the decay ir + — > is: 

2(1 _ p) ( 3 ) 

where r = (m^/m^) 2 . The neutrinos from muon decay fi + — > u fl v e e + have a more implicit 
expression: 

Fv{E v ) = f *!L F,(E,) (g (y) - P^EJ 9l (y)) (4) 

where = .E^/y and gi are known polynomials: go = 5/3 — 3y 2 + 4/3y 3 and 51 = 
1/3 - 3y 2 + 8/3y 3 when v = z^, while y = 2 - 6y 2 + 4y 3 and 3l = -2 + 12y - 18y 2 + 8y 3 
when v = v e The muon flux (from tt + ) that appears in previous formula is: 

F ^ E) ~ 2(1 = ( } 

while the muon polarisation averaged over the pion distribution is given by: 

P (F) x F (F) - gr(gl±jV(gAO , r r ^ 

It is easy to check that in the special case of power law distributions these equations 
reproduce the results of Sect. 7.1 of 0] (e.g., eq. 7.14 there). 

We include the contribution to 7 flux from rj — > 77 (resp., the contribution to v flux 
from the leptonic decay) in the simplest conceivable approximation: namely, we declare 
that the relevant flux of eta mesons (resp., the one of charged kaons) is proportional to 
the one of the neutral pions (resp., of the charged pions) with a fixed coefficient f v = 10 % 

1 Therefore, this procedure is not of direct application for a number of possible galactic sources of 
neutrinos such as micro-quasars |12| that are intrinsically non-transparent or even for extragalactic sources 
since the IR photons background absorbs the VHE gammas above ~ 10 TeV; see also |13|. 
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(resp. fx = 25 % x0.635). Thus: (1) all formulae above should be multiplied by 1/(1 + J^), 
and then (2) we add a neutrino contribution that has the same form as the one from charged 
pions, but replacing r = (m^/mpc) 2 and including the multiplicative factor fx- 
Finally, we incorporate 3 neutrino oscillations replacing: 

~ * Pfifi + F Ve P efJ , (7) 

and the same for antineutrinos. The numerical values of the oscillation probabilities are 
Pnn = 0.39 ±0.05 and P eiM = 0.22 + 0.05 where the quoted errors, approximately equal and 
opposite, are mostly (0.04) due to the spread of 623 around maximal mixing and partly 
(0.02) to the spread of #13 around zero; the effect of the uncertainty in 6*12 is smaller. 
See |14j for more discussion, ^Hj for a resume of neutrino data and analysis, and for 
further references. 

We note in passing a stricter condition on the behaviour of the flux of VHE secon- 
daries with the energy. Consider the connection with the primary cosmic rays F n (E) oc 

dE' F p (E')k(E I E') / E , that we assume for simplicity to be protons. When we go from 
E = E\ to E = E2 with E2 > Ei, the integral decreases because 1) the lower limit in- 
creases; 2) the scaling distribution k in the integrand is a decreasing function; 3) there is 
an explicit factor 1/E. Thus, also f n (E) = EF 7T (E) decreases. The same can be said of 
the function f-^E) = EF 1 {E), since f~/(E) = 2 f n (E/ z)dz; in other words the flux of 
hadronic 7-rays decreases at least as 1/E at high energies. 2 



3 Application: neutrinos from RX J1713. 7-3946 

We apply the formalism of the previous section to obtain the expected neutrino flux from 
RX J1713. 7-3946 on the basis of H.E.S.S. observations [2]. We use 2 parameterisations of 
the 7-ray flux that describe well the observations 2 : 

f 20.4 E- 1 - 98 exp(-£/12) [exponential cutoff] 

F l( E ) ~ j 2Q 1 £_ 2 . 06 ^ + ^ 2 .5]-0.496 [ broken power l aw ] 

where the units are TeV for the energy E and 10~ 12 /(TeV cm 2 s) for the flux F 1 {E). 
We do not use the third parametrization proposed in F~ cx £ , - 2 -°8-o.3iog£;^ name iy 
the distribution with energy dependent photon index: in fact, this cannot result from 
7r°s, since this is just a Gaussian in the logarithmic variable log E, that increases rather 
than decreasing before 40 GeV. Note that a relatively low cutoff implied by H.E.S.S. 
observations is consistent with the present theoretical expectations |18j (however, rather 
different models of the same object are also discussed 19 ) . The result for muon neutrinos 
according to formulae 01 IU and [7| is presented in figureHand tabled i n our approximation, 
the flux of antineutrinos is the same. 



2 Such a very hard spectrum would follow from a hypothetical population of very energetic primaries. 
In fact, consider F P {E') = S(E' - E ): when E < E we find that the pions have Fn(E) oc 1/E since 
fc(0) / 0; thus, the 7— rays would obey the 1/E distribution. 
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Figure 1: Expected Va fluxes corresponding to the two ~f-ray spectra of eq.\$ 
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We can estimate the number of through-going events in a neutrino telescope with 



E th = 50 GeV, <j) = 42°50' (ANTARES location) following ll. Considering neutrinos 
with energies below E umax = 300 TeV (that is not a significant limitation), we find: 3 



that can be compared with N^+p = 8.8 of [TI|, obtained assuming a power law distribution. 
Thus, the new H.E.S.S. data suggest a signal about 8 times weaker than given in [2U] (resp., 
2 times weaker than in ^1]) that adopted a power law distribution extrapolated from the 
first observations of CANGAROO (resp., of H.E.S.S.). 

One can gain something if some events above the horizon are accepted; e.g., with 5° 
more, one can go from a fraction of time useful for observation of 78 % (used for the 
numbers quoted in eq. EJ) to 88 %. This is similar to effects here neglected, e.g., other 
contributions of rj and K meson decays or the deviations of #23 from maximal mixing, 
and should be comparable with the error of the method of calculation we proposed. The 
effect of finite detection efficiency for realistic detector configurations instead should be 
more important (comparable with the effect of the deviation from a power law distribution 
discussed here) for the events are not expected to be particularly energetic: the distribution 
of parent neutrino energies has a median of 3 TeV for both distributions of eq. |H1 

3 The same numbers are obtained using F v + Ft> = 0.37_F 7 . 




4.8 per km 2 per year [exponential cutoff] 
5.4 per km 2 per year [broken power law] 



(9) 
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Table 1: 1 line: selected values of neutrino energy. 2 line: sum of the yields of muons 
and antimuons (including Earth absorption), times the reference area A = 1 km 2 and 
observation time T = 1 year. 3™ d and 4 rd line: the F^^/Fj ratio for the ^y-ray fluxes of 
eq. which varies significantly with the energy. 



E [TeV] 


.1 


.3 


1 


3 


10 


30 


100 


300 


AT(Y^+Y^ [cm 2 s] 


1.0e9 


2.3e10 


2.8E11 


1.9e12 


1.1e13 


3.5e13 


9.2e12 


1.6e14 


F v /F^, exp. cutoff 


.26 


.26 


.25 


.21 


.14 


.06 


.02 


.01 


F v /Ky, brok. power 


.25 


.25 


.25 


.21 


.14 


.13 


.13 


.13 



4 Summary and discussion 

In summary, we presented a simple procedure to convert the observations of high-energy 
7-rays into expectations for high-energy neutrinos, assuming that the source is gamma- 
transparent and that the flux of VHE 7-ray is due to cosmic ray interactions (=it is of 
hadronic origin). The latter hypothesis shows that our flux should be thought as an upper 
bound for gamma-transparent sources. As an application, we calculated the neutrino flux 
from RX J1713. 7-3946 expected on the basis of new H.E.S.S. results and found that the 
expected number of events decreases by 40 — 50 % and that the signal consists of relatively 
low energy events. 

In the future, it will be interesting to repeat the same steps for other intense sources 
of VHE 7-rays, e.g., Vela Jr (RX J0852. 0-4652), that is almost continuously visible from 
ANTARES (95 % of time). In the region E < 10 TeV 1 the spectrum is described by = 
21 E~ 2A (same units as in eq. [SJ. Suppose that the future observations will demonstrate a 
milder exponential cutoff, described by a multiplicative factor exp(— E/E cut ) with E cut = 
50 (150) TeV. In this case we would find Nu+a = 10 (14) events per km 2 per year with a 
median neutrino energy of 5.5 (8.5) TeV (if, again, we assume that all 7-rays are hadronic). 
If instead RX J1713. 7-3946 should turn out to represent a typical SNR in a typical stage, 
it will be important to understand the cosmic rays from a few hundred TeV till the knee, 
e.g., considering other galactic point sources of cosmic rays and/or further phases of cosmic 
ray acceleration. 

These results emphasise even further the importance to obtain 7-ray observations in 
the region from 10 to 100 TeV and to understand well the experimental background coming 
from atmospheric neutrinos. 

We gratefully thank F. Aharonian, V. Berezinsky, P.L. Ghia, D. Grasso and especially 
P. Lipari for useful discussions and M.L. Costantini for help. 
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Note added 

After this work was submitted for publication and when it was presented at Vulcano 
2006 conference (May 2006), a number of interesting new works appeared: where the 
background and possible strategies for neutrino search are quantitatively discussed; [22]) 
where a detailed parameterizations of neutrino and gamma yields is offered; [2H]; where 
neutrinos events from Vela Jr are estimated (though, without describing the details of the 
calculation) using a fixed f/7 conversion coefficient = 1/2. 
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